Draft version December 1, 2011 

Preprint typeset using I^Tf^X style emulateapj v. 08/22/09 



O 

<N 

> 
O 

o 

m 

O 

U 

Oh- 

6 
m 

<N 



THE BIRTH OF A GALAXY: PRIMORDIAL METAL ENRICHMENT AND STELLAR POPULATIONS 

John H. Wise 1 , Matthew J. Turk 2 , Michael L. Norman 3 , Tom Abel 4,5,6 

Draft version December 1, 2011 

ABSTRACT 

By definition, Population III stars are metal-free, and their protostellar collapse is driven by molec- 
ular hydrogen cooling in the gas-phase, leading to large characteristic masses. Population II stars 
with lower characteristic masses form when the star-forming gas reaches a critical metallicity of 
10 -6 — 10 -3,5 Zq. We present an adaptive mesh refinement radiation hydrodynamics simulation 
that follows the transition from Population III to II star formation. The maximum spatial resolution 
of 1 comoving parsec allows for individual molecular clouds to be well-resolved and their stellar as- 
sociations to be studied in detail. We model stellar radiative feedback with adaptive ray tracing. A 
top-heavy initial mass function for the Population III stars is considered, resulting in a plausible distri- 
bution of pair-instability supernovae and associated metal enrichment. We find that the gas fraction 
recovers from 5 percent to nearly the cosmic fraction in halos with merger histories rich in halos above 
1O 7 M0. A single pair-instability supernova is sufficient to enrich the host halo to a metallicity floor of 
1O _3 Z0 and to transition to Population II star formation. This provides a natural explanation for the 
observed floor on damped Lyman alpha (DLA) systems metallicities reported in the literature, which 
is of this order. We find that stellar metallicities do not necessarily trace stellar ages, as mergers of 
halos with established stellar populations can create superpositions of t — Z evolutionary tracks. A 
bimodal metallicity distribution is created after a starburst occurs when the halo can cool efficiently 
through atomic line cooling. 



Subject headings: cosmology 
formation 



methods: numerical — hydrodynamics — radiative transfer — star 



1. MOTIVATION 

The first (Pop III) stars are metal-free and have 
a large characteristic mass. Previously, these stars 
were believed to be massive, with universally sup- 
pressed fragmentation ([Abel et al.l 12002b iBromm et al.l 
l2QQ2HO'Shea fc Norma n 2007). Recently, several studies 
have suggested that while fragmentation at times prior to 
the formation of the first protostar may occur at scales 
of 1000 — 2000 AU, such fragmentation is l ikely to re- 
sult i n high-mass, low- multiplicity groups (jTurk et al.l 
2009; iStacv et "all l2010h . Subsequent studies have sug- 
gested that at later times, subsequent to the formation 
of the first protostar, additional fragmentat ion may pro- 
ceed (jGreif et al.l 1201 lat IClark et al.l 1201 lft . A fraction 
of high- mass, metal- free stars enrich the surrounding 
intergalactic medium (IGM) when they go supernova, 
which can happen in stars < 40 M in Type II super- 
novae (SNe) or in stars roughly betwee n 140 and 260 
in pair-instability SNe (PISNe; iHeger fe Wooslevl 
2002a). The host halo and the neighboring halos are 
then enriched with this eject a. There exists a critical 
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metallicity tha t is ~ 10 6 Z^ if dust coolin g is efficient 
(jOmukai et al.l 120051: iSchneider efall 120061 : IC lark et al.l 



2008) and ~ 10~ 3 5 Z^ otherwise (jBromm et al.l [2001: 



Smit h et al.l l2009). where the gas can cool rapidly, lower- 



ing its Jeans mass. An intermediate characteristic mass 
of ~ 10 Mq can occur if the gas cooling is suppressed 
to the cosmic microwave bac kground (CMB) tempera - 
ture (jLarsonl fl998L iTumlinsonl 120071 : [Smith et al.ll2009l ). 
The resulting Population II star cluster will thus have 
a lower characteristic stellar mass than its metal-free 
progenitors. These first stellar clusters may be con- 
nected to stars in the Milky Way halo and nearby dwarf 
spheroidal (dSph) galaxies, both with a metallicity floor 
of [Z/Hl 7 = -4 dBeers fc Christliebl 120051: iFrebel et al.l 
1201(1 Irafehnever et al.ll2010D . 

Pop III stars primarily form in dark matter halos with 
M ~ 1O 6 M . Their stellar radiation ionizes and heats an 
H II region with a radius 1-3 kpc. The over-pressurized 
H II region drives a 30 km s _1 shock that is 10 times 
greater than the escape velocity of the halo, leaving be- 
hind a warm (3 x 10 4 K) and diffuse (0.5 cm" 3 ) medium 
(iKitavama et al.l 120041 : IWhalen et at] 120041 : lAbel et al.l 
120071 ). This aids in the expansion of the blastwave be- 
cause it delays the transition to the Sedov- Taylor and 
snowplow phases. In PISNe, approximately half of the 
meta ls stay in the IGM with a metal bubbl e size of a few 
kpc flWise fc Abell 120081 : iGreif et al.l l2010l ). The blast- 
wave may induce star formation (SF ) in nearby halos 
through the compression of the gas (|Ferraral fl998). and 
timescales for me tal mixing into the den se gas are many 
dynamical times (|Cen fc R iquelme 2008|) for shock veloc- 



7 We use the conventional notation, 
log(Z /H ). 



[Z/H] = log(Z/H) 
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ities < 100 km s" 1 . 

The transition from Pop III to Pop II SF is 
solely dependent on the propagation of metals from 
the SNe remnants into future sites of SF. Their 
flows are complex because of the interactions be- 
tween the SN blastwave, cosmological accretion, halo 
mergers, and nearby stellar feedback. Spurred 
by the notion of a critical metallicity, this transi- 
tion has been extensively st udied with (!) volume- 
averaged semi-analytic models (fSc annapie co et al.l 2003; 
lYoshida et all \200i IFurlanetto fc Loe b 2005). (ii) mod- 
els us ing hierarchical merg er trees (Tumlinson j 120061 
120101 : ISalvadori et~aTl 120071 : iKomiva et al.l 120101 ). (hi) 
post-processin g of cosmologic al simulations with blast- 
wave models (IKarlsson et al.l 120081 : iTrenti et all 120091 : 
iTrenti fc Shulll l2010f ). and (iv) direct numerica l sim- 
ulations with stellar feedback ([ Tornator e et al.l 120071 : 
iRicotti et all 120081 : iMaio et al.ll2010l l2011f ). 

Over the past decade, these works have refined the 
general picture of inhomogeneous metal enrichment and 
the transition to Pop II star formation, and here we 
give an overview of its development. C onsidering only 
Pop III star formation, Yos hida et al.l calculated that 
Pop III stars can raise the mean metallicity to —4.5 < 
[Z/H] < — 3.5 by redshift 15 in the upper limit where 
all metal-free stars have M = 2OOM and produce 
PISNe. Considerin g both Pop III and II star formation, 
iScannapieco et all found that the transition between the 
two modes is a gradual process where both modes are 
coeval, confirmed by most of the later works. Because 
the host galaxies have small masses and are subject to 
negative radiative feedback through photo-heating and 
photo-dissociation, the minimum halo mass gradually in- 
creases with the radiation background intensity, which 
IFurlanetto fc Loebl found to delay metal enr ichment and 
place it closer to the epoch of reionization. ITrenti et al.l 
noted that underdense regions of the universe are still 
pristine at z = 6 with Pop III stars still forming at these 
late epochs. This group later expanded on these results 
to find that Pop II SFR becomes dominant at z < 25 in 
the buildup of a Milky Way (MW) like halo. Further- 
more they stress the importance of a photo-dissociating 
radiation background that reduces enrichment by PISNe 
and increases the importance of metal-rich galactic out- 
flows, where only 10 -4 — 10 -2 of PISN ejecta ends up in 
extremely metal poor (EMP) stars with M > O.9M . 

The Pop III IMF should play an important role in 
abundance patterns in EMP stars in the MW halo and 
the phys ics of reion ization. Using these data as con- 
straints, iTumlinsonl found that Pop III IMFs with log- 
normal distributions with mean masses between 8 and 
42 M best fit the data. Furthermore he concludes that 
EMP stars with [Z/H] < -3 have between 1 and 10 metal- 
free SN progenitors and the Pop III SFR is less than 1% 
of the total SFR at z = 6. IKarlsson et al.l use observa- 
tional data of EMP stars with their model to constrain 
the mass fraction of Pop III stars that die with a PISN 
is less than 40 per cent. They also conclude that stars 
enriched primarily by PISNe have [Ca/H] > -2.6, which 
could have been missed by some EMP surveys. 

These works have clearly shown that many processes 
play a key role in the transition from Pop III to II stars, 
but how does their interplay affect the formation of the 



first galaxies? Will metal-free pockets of gas still exist in 
dwarf galaxies that are already forming metal-enriched 
stars? How do the metal-enriched ejecta propagate into 
future sites of star formation? Numerical simulations 
are powerful tools to detangle and study these complex- 
ities. To answer such questions, we have carried out 
a simulation that includes both types of SF and their 
radiative and mechanical feedback. The methods used 
here incorporate and link together recent results from 
metal-enriched and metal-free star formation, the criti- 
cal metallicity, and pair- instability supernovae. This is 
the first time it has been possible to correlate the forma- 
tion and feedback of the first stars to protogalaxies, re- 
solving individual molecular clouds and minihalos down 
to 1O 5 M and including the most important physical ef- 
fects, such as radiative and supernova feedback. 

We first describe the simulation setup and numerical 
methods behind our work in ^2) Then in $3] we present 
the evolution of two selected high-redshift dwarf galaxies, 
focusing on their metallicities and stellar populations. 
Last we discuss possible dependencies of our results on 
the Pop III initial mass function and conclude in §4\ 

2. METHOD 
2.1. Simulation setup 

We use the adap tive mesh re finement (AMR) code 
Enzo v2.0 8 ( O'Shea et al.l |2004j ) It uses an TV-body 
adaptive particl e -mesh solver ([E fstathio u et al.l 1985; 
iCouchmanl 119911 : iBrvan fc Normanl Il997h to follow the 
dark matter (DM) dynamics. It solves the hydro- 
dynamical equation using the second-order accurate 
piecewise parabol ic method (| Woodward fc Colellalll98 4l 
IBrvan et al.l ll995). while a Riemann solver ensures accu- 
rate shock capturing with minimal visc osity. We use the 
recently added HLLC Riemann solver (jToro et al.l 1 1994) 
for additional stability in strong shocks and rarefaction 
waves. We use the nine-species (H I, H II, He I, He II, 
He III, e~, H2, H2 + , H ~ ) non- e q uilibr ium chemistry 
model in Enzo (jAbel et al.l 119971: lAnninos et al.l [1 997) 
and the H2 cooling rates from lGlover fc Abell ((2008). We 
spatially distinguish metal enrichment from Pop II and 
Pop III stars. We will follow-up this study with one 
that considers radiati ve cooling fr o m me tal-enriched gas, 
using the method of iSmith et al.l (|2008l ) that is already 
implemented in Enzo. 

To resolve minihalos with at least 100 dark matter 
(DM) particles and follow the formation of the first gen- 
eration of dwarf galaxies, we use a simulation box of 
1 Mpc that has a resolution of 256 3 . This gives us a 
DM mass resolution of 1840 M@. We refine the grid on 
baryon overdensities of 3 x 2 _0,2Z , where / is the AMR 
level, resulting in a supe r-Lagrangian behavior (also see 
Q'Sh ea fc Nor man 2008). We also refine on a DM over- 
density of three and always resolve the local Jeans length 
by at least 4 cells, avoi ding artificial fragme ntation dur- 
ing gaseous collapses ([Truelove et al.l 119971 ) . If any of 
these criteria are met in a single cell, it is flagged for 
further spatial refinement. 

We initialize the simulation with grafic (Bertschinger 
l2001f ) at z — 130 and use the cosmological parame- 
ters from the 7-year WMAP ACDM+SZ+LENS best 

8 enzo.googlecode.com, changeset b86d8ba026d6 
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fit (jKomatsu et all 1201 ID : Vt M = 0.266, ^ A = 0.734, 
Q b = 0.0449, h = 0.71, <r 8 = 0.81, and n = 0.963 with the 
variables having their usual definitions. We use a maxi- 
mum refinement level of / = 12, resulting in a maximal 
comoving resolution of 1 pc. The adaptive particle-mesh 
solver has a force resolution of two cell widths of a given 
AMR grid. We stop the simulation at z = 7 to prevent 
any large-scale modes with r ~ Lbox/2 from entering the 
no n- linear growth regime. At the final redshift, the sim- 
ulation has 1.4 x 10 computational cells and required 
250,000 CPU hours, using 512 cores. 

2.2. Star formation 

We distinguish Pop II and Pop III SF by the total 
metallicity of the densest cell in the molecular cloud. Pop 

II stars are formed if [Z/H] > —4, and Pop III stars are 
formed otherwise. We do not consider Pop III. 2 stars 
and intermediate mass stars from CMB-limited cooling. 

2.2.1. Population III star formation 

Simulations following the full accretion history of Pop 

III stars, from their formation inside cosmological mini- 
halos through transitioning onto the main sequence, are 
not currently possible. However, under the assumption 
of spherical accretion, reasonable estimates of the mass 
ranges for Pop III stars can be made by extrapolating for- 
ward in time the instantaneous accretion rate at the time 
of formation of the first protostellar core and comparing 
the timescales for accretion against the Kelvin-Helmholtz 
time. Several studies, conducting this extrapolation, 
have shown that t he characteristic mass of Pop III stars 
Mchar ~ 100 Ma ( l Abel et al.l 2002; B romm et all 120021 : 
IP 'Shea fe NoTmanl [20071 : lYoshida et all 120081) . They 
form in molecular clouds that coexist with the dark 
matter halo center with limited fragmentation o ccur- 
ring during t heir c ollapse; however iTurk et all ([2009) and 
iStacv etaLl (|201Q[ ) have shown that Pop III binaries may 
form in a fraction of such halos. There have been some 
recent hints that the accretion disk around the first pro- 
tostar may fragment and produce clustered Pop III star 
formation ([Clark et a,l-M2011b IGreif eTa!1l2011al ). The fi- 
nal stellar masses in the aforementioned calculations are 
still highly uncertain with a limited number of samples. 
Nevertheless, they all point toward an IMF with a larger 
characteristic mass than present-day star formation. 

For Pop I I I star s, we use the same SF model as 
IWise fc Abell (|2008f ) where each star particle represents 
a single star. In this model, a star particle forms when a 
cell has all of the following criteria: 

1. an overdensity of 5 x 10 5 (~ 10 3 cm -3 at z = 10), 

2. a converging gas flow (V • v gas < 0), and 

3. a molecular hydrogen fraction /h2 > 5 x 10 -4 . 

These physical conditions are typical of collapsing metal- 
free molecular clouds ~ 10 Myr before the birth 
of a Pop III main-seq uence star (Ab el et all 2002; 

O'Shea fc Normanl 12007 ). This prescription is similar 

to the iCen fc Ostrikern i992) method but removes the 
Jeans unstable (M gas > Mj) and cooling timescale 
(^cooi < ^dyn) requirements. We do not consider the for- 
mer criterion because it is not applicable to simulations 



that resolve the Jeans length at all times. The molecular 
hydrogen fraction requirement effectively constrains star 
formation to cooling molecular clouds, where the H2 for- 
mation rate is significantly larger than the dissociation 
rate from a Lyman- Werner radiation field. 

If multiple cells meet the star particle formation crite- 
ria within 1 pc, we form one Pop III star particle with the 
center of mass of these flagged cells to ensure that one 
massive star is created per metal-free molecular cloud. 
Instead of using a fixed stellar mass like in our previous 
work, we randomly sample from an IMF with a func- 
tional form of 



f (log M)dM = M~ l s exp 



M< 



char 



M 



1.6 



dM (1) 



to determine the stellar mass. Above M c h ar , it 
behaves as a Salpeter IMF but is e xponentially 
cutoff below that mass (Chabriei] l200l IClark et al.l 
2009). For reproducibility, we record the number of 
times the random numbe r generator (Mersenne twister; 
iMatsumoto fc Nishimural ([ 19981 )) has been called for use 
when restarting the simulations. We stress that the Pop 
III IMF is still highly uncertain. In Sec. 2J we discuss 
possible dependencies of our results to changes in the 
IMF. 

After the star particle forms and its mass is randomly 
sampled from the IMF, an equal amount of gas is re- 
moved from the computational grid in a sphere that con- 
tains twice the stellar mass and is centered on the star 
particle. The star particle acquires the mass- weighted 
velocity of the gas contained in this sphere. 

2.2.2. Population II star formation 

We t reat Pop II st a r form ation with the same prescrip- 
tion as IWise fc Cenl ([20091 ), which is similar to the Pop 
III prescription but without the minimum H2 fraction re- 
quirement. This is removed because the metal-enriched 
gas can efficiently cool e ven in the presence of a stron g 
UV radiation field (e.g. Safranek-Shrader et al. 2010). 
To ensure the volume is cooling, we restrict star forma- 
tion to gas with temperatures T < 1000 K. Unlike Pop 
III star particles that represent individual stars, Pop II 
star particles represent a star cluster of some total mass 
and an assumed normal IMF. 

Once a cell meets these criteria, the prescription 
searches outward with increasing radius for the bound- 
ary of the molecular cloud, centered on the most dense 
cell. Here, we define a molecular cloud as a sphere with 
a dynamical time tdyn = 3 Myr (corresponding to an av- 
erage density pMC — 1000/2 cm -3 ) and a radius Rmc, 
where \i is the mean molecular weight. Once the sphere 
radius is found, a fraction c* = 0.07/ co id of the cold gas 
(T < 10 3 K) inside the sphere is converted into a star par- 
ticle with mass m* = c*(47t/3)pmc^mc- This treatment 
of cold gas accretion is similar to a star formation model 
with a multi-phase subgrid model ( Springe ! fc Hernquistl 
l2003f ) but is employed in a simulation that can resolve 
the multi-phase interstellar stellar medium (ISM). 

After the star particle is created, we replace the sphere 
with a uniform density pMC = (1 — c *)/(^d yn ) an d 
temperature T = 10 4 K, which approximates the initial 
stages of an H II region. Our choice of c* agrees with the 
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theory proposed in iKrumholz fc McKeel (|20Q5f ) that are 
supported by observations. Observations of star- forming 
regions show that the dynamical time is approximately 
700 kyr an d that SF is ex tended over several dynamical 
times (e.g. iTan et aLll2006f ). which motivates our choice 
of tdyn of our definition of molecular clouds in our SF 
prescription. 

Similar to the Pop III treatment, we merge any newly 
created particles within a radius of 10 pc. The new star 
particle acquires the center of mass of the merged par- 
ticles. During the merging, the gas is unaltered besides 
the initial gas accretion from the grid to the star par- 
ticles. We set the minimum mass of a star particle to 
^*,min = 1OOOM . If the initial mass does not exceed 
^*,min, the star particle does not provide any feedback 
and continues to accrete until it reaches ra* 5m i n . 

2.3. Stellar feedback 

The radiation field is evolved with adapti ve ray tracing 
(|Abel fc Wandemi200alWise fc Abelj[20lTI) that is based 
on the HEALPix framework ([Gorski et al.ll2005f ) and is 
coupled self-consistently to the hydrodynamics. In this 
work, we only consider energy coupling to the gas; we 
defer an investigation of the effects of radiation momen- 
tum coupling to a subsequent paper. Each star parti- 
cle is a point source of hydrogen ionizing radiation with 
the ionizing luminosity equally split between 48 initial 
rays (HEALPix level 2). We use a mono-chromatic spec- 
trum for the radiation with the energy equaling the 
luminosity- weighted photon energy of the spectrum. For 
a cosmo logical simulation that focuses on galaxies, this 
does not signifi cantly affect the ove rall galactic dynamics 
(see Sec. 6.3 in I Wise fc Abell 12011). We do not consider 
helium singly- and doubly-ionizing radiation. As the rays 
propagate from the source or into a high resolution AMR 
grid, they are adaptively split into 4 child rays, increas- 
ing the angular resolution of the solution, when the solid 
angle of the ray = 47r/(12 x 4 l ) is larger than 1/3 of the 
cell area, where L is the HEALPix level. We model the 
H2 dissociating radiation with an optically-thin, inverse 
square profile, centered on all Pop II and III star parti- 
cles. In principle, self-shielding effects can be included by 
utilizing our ray tracing framework or other approximate 
methods that calculat e the H2 column den sity using the 
local Sobolev length (Wol cott- Green et al .1120 ill ) or only 
i ntegrating along Cartesian axes or random lines of sight 
(lYoshida et al. 1 120031 [20071 : 1 Glover fc Mac Lowll2007l) . In 
halos with M > 1O 8 M , the H2 column density may be- 
come large enough t o self-shield itself from LW radiation 
(jWise fc A bel 2008). Thus we may be underestimating 
the SFRs in our work. 

2.3.1. Population III stellar feedback 

We use mass-dependent hydrogen ionizing and LW 
phot on luminosities and lifetimes of the Pop III stars 
from |Schaerer (2002). We use a mass-independent pho- 
ton energy E p h = 29.6 eV, appropriate for the near- 
constant 10 5 K surface temperatures of Pop III stars. 
They die as Type II SNe i f 11 < M*/M < 40 
M (IWooslev fc Weaverl[l995h and as PISNe if t hey are 
in the mass range 140-260 M (jHeger et al.ll2QQ3h . where 
M* is the stellar mass. For normal Type II SNe be- 
tween 11-20 Mq, we use an explosion energy of 10 51 erg 



and a linear fit to t he metal eject a mass calculated in 
iNomoto et~aT1 (|2C 



M z /M = 0.1077 + 0.3383 x (M*/M - 11). (2) 

We model the SNe of stars with 20 < M*/M < 40 
as hyperno ya with the en ergies and ejecta masses also 
taken from lNomoto et all linearly interpolating their re- 
sults to_J\| i _ : _J^o^^ the explosion energy 
from He ger fc Woosleyl (j2002bf ). where we fit the follow- 
ing function to their models, 

£pisn = 10 51 x [5.0 + 1.3O4(M He /M - 64)] erg, (3) 

where M He = (13/24) x (M* - 2O)M is the helium 
core (and equivalently the metal ejecta) mass and M* 
is the stellar mass. If the stellar mass is outside of these 
ranges, then an inert, collisionless black hole (BH) par- 
ticle is created. The blastwave is modeled by injecting 
the explosion energy and ejecta mass into a sphere of 10 
pc, s moothed at its su rface to improve numerical stabil- 
ity (Wise & Abel 120081 ). Because we resolve the blast- 
wave relatively well with several cells across at its ini- 
tialization, the thermal energy is converted into kinetic 
energy and agre es with the Sedov- Taylor solution (e.g. 
IGreif et alJl2QQ7h . 

2.3.2. Population II star formation 

The Pop II stars emit 6000 hydrogen ionizing photons 
per baryon over their lifetime and i£ p h = 21.6 eV, ap- 
propriate for a [Z/H] = —1.3 population. We note that 
lower-metallicity stars generate up to 60% more ioniz- 
ing photo ns and might b e underestimating the radiative 
feedback (jSchaererl l2QQ3h . Nonetheless this study pro- 
vides an excellent first insight in the transition to Pop II 
SF, as the metal enrichment is the key ingredient here. 
The star particles live for 20 Myr, the maximum lifetime 
of an OB star. These stars generate the majority of the 
ionizing radiation and SNe feedback in stellar clusters, 
thus we ignore any feedback from lower mass stars. 

For SN feedback, these stars generate 6.8 x 10 48 erg 
s _1 M _1 after living for 4 Myr, which is injected into 
spheres with a radius of 10 pc. However if the resolution 
of the grid is less than 10/3 pc, we distribute the energy 
into a 3 3 cube surrounding the star particle. In the same 
region, star particles also return ejected material with 



Arn e 



0.25 At M* 
t - 4 Myr 



(4) 



back to the grid at every timestep on the finest AMR 
level. Here to = 20 Myr is the lifetime of the star par- 
ticle. This ejected gas has solar metallicity Z = 0.02, 
resulting in a total metal yield y = 0.005. In addition 
to numerical stability, energy distribution across several 
cells instead of a single cell has been shown to match cos- 
mic SFRs better in galaxy simulations at lower redshifts 
(jSmith et al.ll2Qllh . 

3. RESULTS 

Here we first present the global star formation rate 
(SFR) during early galaxy formation and then focus on 
the gaseous and stellar evolution of two selected halos in 
the simulation: one that has an early mass buildup but 
no major mergers after z = 12, and one that experiences 
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Fig. 1. — Evolution of the entire simulation volume (Lbox — 1 Mpc) at redshifts 15, 12, 10, 8, and 7 (left to right). Pictured here are 
the density- weighted projections of density (top), temperature (middle), and metallicity (bottom). Note how the stellar radiative feedback 
from low-mass galaxies reionize the majority of the volume. The metallicity projections are a composite image of metals originating from 
Pop II (red) and III (blue) stars with magenta indicating a mixture of the two. 
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TABLE 1 

Enrichment Fractions at Redshift 7 



1Q -27 1Q -2 

Density [cm" 



10 3 10 4 
Temperature [K] 



[Z/H] 


/.(Pop III) 


/ m (Pop III) 


/.(Pop II) 


/m(Pop II) 


-6 


3.81(-2) 


1.43(-1) 


5.52(-3) 


5.22(-2) 


-5 


3.28(-2) 


1.24(-1) 


3.92(-3) 


4.36(-2) 


-4 


2.65(-2) 


1.01(-1) 


2.33(-3) 


3.43(-2) 


-3 


1.86(-2) 


6.45(-2) 


1.00(-3) 


2.50(-2) 


-2 


8.29(-3) 


1.23(-2) 


2.61(-4) 


1.73(-2) 



-6 -4 -2 
[Z 2 /H] 

Fig. 2. — Density- weighted projections of density (left), temper- 
ature (center), and metallicity (right) of the two selected halos at 
z = 7. The field of view is 5 proper kpc, and the circles have a 
radius of r2oo- The composite metallicity images are constructed 
in the same fashion as Fig. [T] 



a series of major mergers between z = 10 and z = 7. We 
name the halos "quiet" and "intense" , respectively. The 
entire simulation contains 38 galaxies with 3640 Pop II 
stellar clusters and captures the formation of 333 Pop III 
stars. 

We illustrate the evolution of the simulation from 



Note. — f v and f m are volume filling and mass- weighted frac- 
tions, respectively, above several different metallicity cuts. Values 
are listed in scientific notation with the exponents in parentheses. 



z = 15toz = 7in Figure [T] with density weighted pro- 
jections of gas density, temperature, and metallicity of 
the entire box. Figure [2] shows the same types of plots 
but focusing on the two halos of interest with a field of 
view of 5 proper kpc. Radiative and mechanical feedback 
create a multi-phase medium inside these halos, which 
are embedded in a warm and ionized IGM. At the final 
redshift, the volume and mass ionized fraction are 0.76 
and 0.72, respectively, which are consistent wi th large- 
scale reionization simulations (Za hn et al.ll201l| ). PISNe 
enrich 0.065 and 0.019 of the mass and volume, respec- 
tively, above 10 -3 Z , and 0.025 and 0.0010 of the mass 
and volume are enriched by Type II SNe from Pop II stars 
above the same threshold. With M c h a r = 100M Q , a neg- 
ligible fraction of Pop III stars die as Type II SNe; thus, 
all metals originating from Pop III stars are produced in 
PISNe. Table [1] gives additional enriched fractions with 
different metallicity cuts. 

3.1. Global star formation rate 
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Time [Myr] 

Fig. 3. — Cosmic star formation rates for Pop II (solid) and III 
(dashed). Note that the Pop III rates are nearly constant from 
redshift 12 to 7. Pristine halos still survive in underdense regions, 
supporting late-time metal-free star formation. 

Figure [3] shows the cosmic SFR in the simulation for 
Pop II and III stars. Metal- free stars form at a nearly 
constant rate of 5 x 10 -5 M yr _1 Mpc -3 until the end of 
the simulation. The dissociating and ionization radiation 
from nearby Pop III stars and galaxies do not halt the 
formation of such stars in low-mass halos but only sup- 
press them. We expect a uniform H2 di ssociating radia- 
tion background to further delay them (Macha cek et al.l 
120011 : iWise fe AbelH2007t IQ'Shea fe Normanl[2QQ8l) ~ 

The cosmic Pop II SFR is best fit with an exponential 
exp(— z/l.l). Metal-enriched star formation begins early 
in the simulation at z ~ 16 when a PISN blastwave over- 
takes a satellite halo with M ~ 5 x 1O 5 M , triggering 
star formation there. At z > 10, these enriched miniha- 
los cause the Pop II SFR to be higher than this relation. 
Only after z ~ 10 does the SFR follows this exponential, 
which is a steepe r time dependency than the exp(— z/S) 
relation found in iHernquist fc Springell (|2003f ) , who also 
acknowledged that this behavior should occur at early 
times. At redshift 7, it reaches 2 x 10 -2 M yr _1 
Mpc -3 . This value agrees with derived SFRs from z = 
6 — 1 galaxy luminosity functions (e.g. iBouwens et al.l 
120111) and large-sc ale galaxy formation simulations (e.g. 
iSchave et al"1l2010l ). 

3.2. Mass evolution 

Figure H^i shows the total, metal-enriched stellar, and 
gas mass history of the most massive progenitors of both 
halos. The quiet halo undergoes a series of major mergers 
at z > 12, growing by a factor of 30 to 2.5 x 1O 7 M 
within 150 Myr. Subsequently, it only grows by a factor 
of 3 by z = 7 mainly through smooth accretion from 
the filaments and IGM. It is the most massive halo in 
the simulation between redshift s 13 and 10. At the same 
time, the intense halo has a mass M = 3 x 1O 6 M , but 
it is contained in a biased region on a comoving scale of 
50 kpc with - 25 halos with M ~ 1O 6 M . After z = 10, 
these halos hierarchically merge to form a 1O 9 M halo 
at z = 7 with two major mergers at redshifts 10 and 
7.9, seen in the rapid increases in total mass in Figure 
HJl. The merger history of the two halos are not atypical 
as dark matter halos can experience both quiescent and 



vigorous mass accretion rates. 

Both halos start forming Pop II stars when M = 
1O 7 M . This is consistent with the filtering mass 
Mj, which is plotted as a thick line in Figure 2^, 
of high-redshift halo s when it accrete s mainly from 
a pre-heated IGM dGn edin fc Huil 11995 IGnedinl [200(3 : 
IWise fc Abell 120081 ). Afterwards these halos can cool 
efficiently through H2 cooling, sustaining constant and 
sometimes bursting SF. The latter characteristics are 
equated with the definition of a galaxy. The quiet halo 
forms 1O 5 M of stars by z = 9. This initial starburst 
photo-evaporates the majority of its molecular clouds, in 
addition to heating and ionizing the surrounding IGM 
out to a radius of 10-15 kpc at z = 9. These respectively 
reduce the in-situ and external cold gas supply that could 
feed future SF. 

The gas fractions of both halos decrease from 0.15 to 
0.08 by outflows driven by ionization fronts and blast- 
waves in their initial starbursts. The quiet halo does not 
have a major merger with any halo with M > Mf, lead- 
ing to a small final gas fraction. These low-mass halos are 
photo-evaporated, hosting diffuse warm gas reservoirs in- 
stead of cold dense cores. After z = 10, the halo mainly 
accretes warm diffuse gas from the filaments and IGM. 
In contrast, the intense halo grows from major mergers 
of halos with M > Mf . The progenitor halos involved in 
the major mergers are able to host molecular clouds and 
have higher gas fractions. Between z = 10 and z = 8, 
the gas fraction increases from 0.07 to 0.12 until it jumps 
to 0.14 when a gas-rich major merger occurs. The stellar 
mass accordingly increases with the ample supply of gas 
during this period. 

3.3. Metallicity evolution 

The evolution of the stellar and gas metallicity of both 
halos are illustrated in Figure HJd. PISNe from Pop III 
stars enrich the nearby IGM out to a radius of 10 kpc 
and provides a metallicity floor of [Z3/H] ~ —3. Ejecta 
from Pop II SNe initially enrich the ISM of both halos 
to an average [Z2/H] between -3.5 and -3. 

In the quiet halo, an equilibrium of [Z2/H] ~ —2.5 is 
established between metal-rich outflows and metal-poor 
inflows. Galactic outflows are directed in the polar di- 
rections of the gas disk, keeping the adjacent filaments 
metal-poor. These features and a well- mixed ISM (cf. 
IWise fc Abeil2QQa iGreif et al.ll2010l ) are apparent in the 
metallicity projections in Figure_[2] The average stellar 
metallicity is within 0.5 dex of [Z2/H]. 

In the intense halo, the first few Pop II star clusters 
have [Z/H] between -1 and -2 and dominate the aver- 
age stellar metallicity at z > 8._ The metallicity later 
increases by a factor of 30 to [Z2/H] = —1.5 through 
self-enrichment from a starburst. Because this halo is 
located in a large-scale overdensity, most of the ejecta 
falls back into the halo after reaching distances up to 
20 comoving kpc, keeping the halo gas metallicity high 
because the inflows are relatively metal-rich themselves. 
After the z ~ 8 starburst, the average stellar metallicity 
follows the average gas metallicity within 0.1 dex. 

To further demonstrate the metallicity floor in high- 
redshift dwarf galaxies, we show the average metallicity 
from Pop III SNe [Z3/H] in all halos as a function of 
halo mass at z = 7, 8, 9, and 10 in Fig. [5j and it is 
strongly dependent on halo mass. The [Z3/H] error bars 
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Fig. 4. — (a) Evolution of the total halo mass (top), stellar mass (middle), and gas fraction (bottom) of the quiet (dashed) and intense 
(solid) halos. In the top panel, the filtering halo mass is plotted as the thick black line (overlapping with the quiet halo after z = 9. (b) 
Mass- weighted stellar metallicities and gas metallicities enriched by Pop II and Pop III SNe of the intense (top) and quiet (bottom) halos. 
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Fig. 5. — Average halo metallicity from Pop III SNe as a function 
of total halo mass. Inhomogeneous metal enrichment causes the 
large spread in metallicity at low masses. At M^aio > 3 x 10 7 Mq, 
all halos are enriched by Pop III SNe above [Z/H] > -4. The y- 
axis error bars show the 25% and 75% percentile values because 
the underlying distribution is not Gaussian. 

mark the 25% and 75% percentiles of the distribution 
at M = Mhaio because the underlying distribution is 
not necessarily Gaussian at all halo masses. Above 3 x 
10 7 M o , all halos enriched to [Z 3 /H] > -4 with little 
time evolution, further illustrating the metallicity floor. 
Below this halo mass, the distribution is roughly bimodal 
composed of pristine and metal-enriched halos. 

3.4. Star formation history 

The most massive progenitor of the quiet halo inter- 
estingly never hosts a Pop III star. Instead a nearby 
halo forms a Pop III star, which (randomly) produces a 
PISNe at z = 16. The blastwave overruns the most mas- 
sive progenitor, and the dense core survives this event 
and is enriched by this PISN, triggering the transition 
to Pop II SF. Other progenitors host three Pop III stars, 
forming at z = 15.4, 14.2, 13.8, with the latter producing 
a PISN. Metal enrichment from these two PISNe and 
Pop II SF quench Pop III SF in this halo. The progen- 



itors of the intense halo host a total of 56 Pop III stars 
with 21 producing PISNe and the other 35 producing 
intermediate-mass BHs. The first forms in a 6 x 1O 5 M 
halo at z = 19. Pop III stars form on a regular inter- 
val in the halo's progenitors until z = 9 when most of 
these halos enter the metal-rich bubble surrounding the 
intense halo. 

Figure [6] shows the SF history (SFH), metallicity distri- 
bution, and SF rates of both halos. A nearby PISNe pro- 
vides a metallicity floor of [Z/H] = -2.8 in the quiet halo 
at which metallicity the first Pop II stars form. The stel- 
lar metallicity evolution exhibits what is expected from 
an isolated system with the stellar feedback steadily en- 
riching the ISM, resulting in a correlation between stel- 
lar age and metallicity. After z = 10 the metallicities 
plateau at [Z/H] = -2.1 for reasons previously discussed. 
The SFR peaks at z = 10 and decreases as the cold 
gas reservoir is depleted. Around z = 7.5, a 25:1 mi- 
nor merger occurs, and the gas inside the satellite halo 
is compressed, triggering metal-poor, [Z/H] = -3.2, SF 
during its nearest approach. This halo remains metal- 
poor because most of metal enrichment in the quiet halo 
occurs in bi-polar flows perpendicular to the galaxy disk 
and filament. Stars with [Z/H] < -3 compose 1.6 percent 
of the total stellar mass. 

In contrast with the quiet halo, the intense halo un- 
dergoes a few mergers of halos with an established stel- 
lar population. This creates a superposition of age- 
metallicity tracks in the SFH, seen in the complexity 
of Figure [6l The first two Pop II stellar clusters have 
an unexpectedly high metallicity [Z/H] ~ -1, which oc- 
curs when a PISN blastwave triggers SF in two neighbor- 
ing halos. Most of the early SF have [Z/H] = -2.5. At 
z = 9, the halo's virial temperature reaches 10 4 K. This 
combined with a 10:1 merger creates a starburst that 
quickly enriches the halo to [Z/H] = —1.5 by z = 8. The 
halo continues to enrich itself afterwards. The spikes in 
the scatter plot correspond to SN triggered SF in nearby 
molecular clouds that are enriched up to a factor of 10 
with respect to the ISM. However their mass fraction are 
small compared to the total stellar mass. The starburst 
at z = 9 creates a bimodal metallicity distribution with 
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Fig. 6. — The scatter plots show the metal-enriched (Pop II) star formation history of the intense (left) and quiet (right) halos as a 
function of total metallicity, i.e. the sum of metal ejecta from both Pop II and Pop III SNe, at z = 7. Each circle represents a star cluster, 
whose area is proportional to its mass. The open circles in the upper right represent 10 3 and 10 4 Mq star clusters. The upper histogram 
shows the SFR. The right histogram depicts the stellar metallicity distribution. The intense halo shows a large spread in metallicity at 
z > 10 because these stars formed in progenitor halos that were enriched by different SN explosions. At z < 10, the majority of stellar 
metallicities increase as the halo is self-enriched. The spikes in metallicity at t = 620, 650, and 700 Myr show induced star formation with 
enhanced metallicities in SN remnant shells. The dashed lines in the left panel guide the eye to two stellar populations that were formed 
in two satellite halos, merging at z = 7.5. The quiet halo evolves in relative isolation and steadily increases its metallicity to [Z/H] ~ -2 
until there is an equilibrium between in- situ star formation and metal-poor inflows from filaments. 




log M tot [M s ] 

Fig. 7. — Total Pop II luminosities of halos versus total halo 
mass with constant mass-to-light ratios overplotted at redshift 7. 
The halos with M > 1O 8 M0 form stars efficiently because they can 
cool through atomic line cooling, whereas the lower mass galaxies 
rely on molecular hydrogen cooling only. The low M/L ratios and 
spread across two orders of magnitude are similar to ones observed 
in nearby ultra faint dwarf galaxies. 



peaks at [Z/H] = -2.4 and -1.2 with the metal-rich com- 
ponent mainly being created after the starburst. Two 
systems with sizable stellar components merge into the 
halo at z ~ 8, and their stellar populations, which are 
traced by dashed lines in Figure El are still discernible in 
the metallicity-age plot. Stars with [Z/H] < -3 compose 
1.8 percent of the total stellar mass. 

3.5. Mass-to-light ratios 



Mass-to-light ratios are crucial in interpreting obser- 
vations when connecting the galaxy to its DM halo, and 
ultimately to models of galaxy formation. For this paper, 
we use a simple model to calculate the total stellar lumi- 
nosity of each galaxy. In a later paper, we will follow up 
with more sophisticated models to determine more ac- 
curate observables. Stars in each particle are treated as 
an instantaneous burst. Using the empirical law of main 
sequence lifetime tyis M -5 / 2 , the maximum stellar 
mass in such a system is dependent only on its age so 
= (t*/10 Gyr) -2 / 5 . Then combining another scal- 
ing law L oc M 7 / 2 and the Salpeter IMF ip(M) oc M~ a , 

pM up 

L/Mq — \ i/;(M)L(M)dM 

JM lo 

^M u 9 p / 2 -«(9/2-a)-\ (5) 

where we set M\ Q = because the total luminosity of 
low-mass stars is negligible. Then the galaxy luminosity 
is the sum of all star particle luminosities, which we plot 
in Fig. [3 with respect to halo mass, along with the DM 
halo mass function. 

Recall that the most massive halo (M to t = 1O 9 M ) is 
undergoing a starburst at z = 7, translating into a rel- 
atively high M/L = 6. The lower mass galaxies show a 
spread in M/L between 30 and 3000, similar to the Sloan 
Digital Sky Survey (SDSS) ob servations of local ultra 
faint dwarf (uFd) galaxies (e.g. iStrigari et al1l2008f ). It 
has been suggested that these low- mass (M tot < 1O 8 M ) 
galaxies that cannot cool through hy drogen line cool- 
ing are the progenitors of uFd galaxies (jBovill fc Ricottil 
l2011a ,b). where all of the star formation occurs before 
reionization and was subsequently terminated by UV 
heating from external sources. These galaxies have stel- 
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lar masses ranging from 210 M© (SDSS J1058+2843) to 
3.2 x 1Q 5 Me?) (Dra| with most between 10 3 and 10 4 M . 
(jMartin et al.ll2008h . Because the total stellar population 
would be sampled, on average, by several star particles, 
we cannot make robust conclusions about the nature of 
uFd galaxies. 

4. DISCUSSION AND CONCLUSIONS 

We have focused on the birth of two galaxies prior to 
reionization with a cosmological AMR radiation hydro- 
dynamics simulation with very high spatial (1 pc comov- 
ing) and mass resolution (1800 M©). With this resolu- 
tion, we are able to follow the Pop III star formation 
and radiative and chemical feedback effects in hundreds 
of minihalos that contribute to the formation of a proto- 
galaxy. Supernovae from Pop III stars provide the nec- 
essary heavy elements for the transition to a Population 

II stellar population, which we have directly simulated. 
The galaxies simulated here represent a population of un- 
characterized high-redshift dwarf galaxies, which should 
be detected with the upcoming James Webb Space Tele- 
scope, that should contribute a significant fraction of 
the ionizing photons during reionization. Recently, a 
z = 8.6 spectroscopically confirmed galaxy was found 
to reside in an H II region that extends beyond 1 Mpc. 
However, using the SFRs and Lya continuum as esti- 
mators, this galaxy produces enough ionizing radiation 
to create an H II region between - 0.1(/ esc /0.1) 1/3 Mpc 
and 0.5(/esc/0-l) 1 ^ 3 Mpc, where f esc is t he fraction of 
ionizing radiation escaping into the IGM (Leh nert et al.l 
l2010f ). This discrepancy could be caused by undetected 
satellite dwarf galaxies, such ones presented here, that 
contribute to local radiation field, or an ionizing ra- 
diation escape fraction near unity ([Wise & Cen 2009; 
i Razoumov fc Sommer-Larsenll20fol : iPaardekooper et al.l 
120111 : lYaiima et al.ll2011h ~ 

We find that one PISN is sufficient to enrich the star- 
forming halo and surrounding ~ 5 kpc to a metallicity 
of 1O _3 Z , given M c h ar = 100 M . DLA systems have 
a metallicity floor on the same order (jWolfe et al.l 12005: 
iPenprase et al.ll2010f h and metal enric hment from Pop 

III SN e provides a possible explanation (jKobavashi et al.l 
l2011f ). Our simulation strengthens this claim of a metal- 
licity floor from Pop III SNe (see also iTornatore et al.l 
2007; Karls son et al.l 120081 : iMaio et al.ll20TTf T where our 
work improves on previous studies by resolving star- 
forming molecular clouds and employing more realistic 
feedback physics, in particular, radiation transport. 

If the first stars have a lower cha racteristic mass that 
favor hypernovae (jTumlinsonl l2Q07f ) . then this metallic- 
ity floor should be lowered by a factor of ~ 10 be- 
cause th^_nieJaJ_ejectafrom a 40 M hyper nova is 8.6 
Mq ( Nomoto et a l. 2006), co mpared to 85 of metal s 
produced by a 180 M© PISN (jHeger fc Wooslevl 12002 of ). 
Both explosions have approximately 3 x 10 52 erg. There- 
fore the blastwaves from these two explosions should ex- 
pand to equal radii before stalling in the IGM. We note 
that stronger metal-line cooling in the PISN case may 
cause the blastwave to stall at a marginally smaller ra- 
dius. 

Building upon this idea of a lower metallicity floor, 
consider a pre-galactic halo that has hosted a Pop III 
star, producing a SN. The majority of the gas is ejected 



i nto the IGM and re- accretes after y 10 — 50 Myr 
(j Johnson et al.ll2007l : IWise fc Abelll2008f ). What happens 
if the metallicity is not greater than the critical metal- 
licity? Massive Pop III stars would form again in this 
halo; however, enriching the surrounding IGM with Pop 
III SNe and ejecting mat erial becomes more difficult in 
halos with M > 1O 7 M (jWhalen et all 12008( 1. Thus it is 
more likely for these halos to retain more of its gas dur- 
ing star formation episodes and self-enrich its host halo 
and less of the IGM. This would accelerate the transition 
from Pop III to II in more massive halos. 

We conclude that it only takes one, at most two, SNe 
from Pop III stars in the halo progenitors to complet e 
the transition to Population II (Frebe l fc Brom m 2010). 
However, this does not prohibit more than two instances 
of Pop III star formation in halo from happening be- 
cause the stars can produce BHs with little or no metal 
ejecta, keeping the surrounding environment pristine. 
Even though we find that Pop III SNe create a floor 
of ~ 10 _3 Z©, gas with lower metallicities can still ex- 
ists in nearby halos, where the metal-rich SN ejecta will 
mix sl owly into dense cores as the blastwave overtakes 
them (jCen fc Riquelm e 2008). These environments may 
be the host of true extremely metal-poor star formation, 
such as the recently disc overed halo star wit h all detected 
metals below 10~ 42 Z© (jCaffau et al.ll2011f ). 

We also find that the merger history plays an im- 
portant role in supplying gas into halos after its first 
epoch of SF. Mergers of halos below the filtering mass 
are inefficient in providing gas whereas the opposite is 
true for merging halos above the filtering mass. Ha- 
los do not necessarily need T v i r > 10 4 K to form sig- 
nificant stellar populations; however SFRs dramatically 
increase, and thus metal enrichment, when this thresh- 
old is reached. Halos with mostly gas-poor mergers or a 
quiescent merger history result in a monotonic increase 
in metallicity with time. SFHs become more complex 
with multiple metallicity-age tracks if the halo experi- 
ences mergers with other halos that have an established 
stellar population. Furthermore SF that is triggered by 
blastwaves interacting with molecular clouds can have 
metallicities up to a factor of 10 higher than the main 
starburst. Starbursts at T v i r = 10 4 K can enrich the host 
halo enough to create a bimodal metallicity distribution, 
where the metal-poor component is created before the 
burst. Note that mergers of stellar populations can also 
create a si milar bimodal d istribution, which are observed 
in dSphs ([Battaglia et al.ll201lD . 

Because our simulation only samples 1 comoving Mpc 3 , 
it misses the more rare densit y fluctuations and is b iased 
toward late galaxy formation (iBarkana fc L oeb 2004J). In 
such rare galaxies, the physical scenario described by 
this work is largely unchanged but is shifted to higher 
redshifts. Thus our results can be generalized to any 
epoch prior to cosmic reionization. However, two possi- 
ble consequences in very high (z > 20) redsh ift galaxy for- 
mation are 1) an intermediate mas s IMF (|L arson 1998; 
iTumlinsonl 120071 : iSmith et al1l2009[ ). resulting from the 
CMB temperature floor, could be more prevalent, and 2) 
stream ing; gas velocities arising from the recombination 
epoch (Tseli akhovich fc Hiratall20Tol) could suppress star 
formation in halos belo w 10 6 M© (jTseliakhovich et al.l 
[2QlTl IGreif et alJl20lTbh . 
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We have shown that it is possible to simulate the for- 
mation of a high-redshift dwarf galaxy and its entire SFH 
with radiative and mechanical feedback. These results 
provide invaluable insight on the first galaxies and the 
role of metal- free stars in the early universe. We plan to 
follow up this work by focusing on observational connec- 
tions with high-redshift galaxies and Local Group dwarf 
galaxies in the near future that can better constrain the 
physical models and assumptions, such as the critical 
metallicity, radiation backgrounds, and IMFs, used in ab 
initio dwarf galaxy simulations. 
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